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ON CERTAIN INTEGRALS WHICH PERTAIN TO THE FORCED VIBRATION OF PLATES 


By 

Harold Levine 

Joint Institute for Aeronautics and Acoustics 
Stanford University 


Introduction 

The normal velocity anywhere on an infinite plane elastic surface 
immersed in a uniform compressible fluid and driven by a time-periodic 
concentrated force or moment admits a single complex Fourier type integral 
representation; and its estimate at points remote from the site of excitation 
can be achieved with the standard procedures of asymptotic analysis. To 
investigate the surface response near the driving point, Crighton (1972) 
employs a noteworthy albeit indirect procedure which (in the case of a line 
force or moment) rests upon a correlation between asymptotic behaviors of the 
respective functions 

CO 

f(x) = i d? , |x| 0 (1) 

<-.00 

and 

00 

F (?) = f f(x)e^^^ dx , |?| -> 00 , Im ? > 0 (2) 

0 

in the designated limits. A particular consequence of the latter relation, 
namely 

00 

F (in) = / f(x)e dx , n > 0 (3) 

0 

associates (through Watson's lemma) terms in the asymptotic expansion of F_^(in) 
as n +°° with the Laplace transform of terms in the expansion of f (x) as 
X -> +0 . The determination of F_j_(?) , or half-range Fourier transform of 
f(x) , is accomplished by a split of the full range transform 
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F(0 = F^(C) + F_(C) 


(4) 


into component parts which are analytic functions of the complex variable C 
in the overlapping domains Im ? > -e , < £ , respectively, F(C) itself being 
analytic within the strip |lm c| < £ of the ^-plane. As Crighton observes, 
however, to effect the above split (which is a characteristic of the Wiener- 
Hopf technique for solving integral equations) generally poses an integration 
problem with comparable difficulty to that presented at the outset by (1) ; 
and there are additional technicalities requiring consideration if (1) is re- 
placed by the Fourier-Bessel integral (or Hankel transform) 


f (x) 


f F(c)J.(?x)C dC 

0 


(5) 


that occurs in the case of point excitation. Prior to Crighton’ s suggested 
function theory approach and with a wholly different motivation, Grosjean 
(1965) gave precise rules and conditions for the explicit construction of 
convergent infinite series, around x = 0 , which describe functions defined 
by the Fourier cosine/sine integrals. 


/F(5) cos 5x d? and /F(^) sin d^ ; 

0 0 

since these are predicated on the hypothesis that F(^) is a real function, 
with a convergent or asymptotic expansion of the form 


a- a« 0 < a ^ 1 

^oH-1 ^ » 

C C a real 

^ n 


the results do not serve the purposes of a near field plate excitation analysis. 


It is the objective of this note to directly recast the various plate 
integrals considered by Crighton in a manner that permits their systematic 
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development near the source point; and thus to dispense with manipulation of 
the velocity transform. The advantage of such a direct procedure is manifest 
by the ease with which distributed forces over the plate are incorporated into 
the analysis, whereas Crighton’s approach is rendered more difficult at the out- 
set by any change in the velocity transform or the function F(?) . After the 
procedure is introduced (§1) in connection with a line source problem and followed 
(§2) by consideration of a point source, the analogous problems relative to uni- 
formly distributed forces over strips or circular areas of the plate receive 
attention (§3). There are good reasons, since the exact plate admittances pos- 
sess a complicated analytical nature, to seek out useful (simple) approximate 
representations along the lines taken by Crighton, in fact, for the heavy fluid 
loading (low frequency) limit with a point force or moment. 


1. The Line Source Problem 

The surface velocity v(x)e of a thin elastic plate lying in the 

(x,z)-plane, with a uniform compressible fluid adjacent thereto (in y > 0 , 
say), which results from the application of a normally directed line force 
along the z-axis, is represented by the Fourier integral 


v(x) 


icjF, 


2ttB 


00 


f 

mmOO 


t^-k^ exp(-l^x) dg 
P P 


( 6 ) 


wherein m and B denote the specific mass and bending stiffness of the 

plate, p = ym is the fluid density, k = 0)/c specifies the acoustic wave 

2 1/4 

number in the fluid, and k^ = (mo3 /B) designates the bending wave number 
of the plate in vacuum. Following a customary prescription, the path or con- 
tour of the integral (6) approaches the real axis as |?| and passes 

above/below all singularities on the negative/positive real axis; and the 

/2 2 

function /? -k , rendered single-valued with cuts from k to <» and -<» 
to -k , respectively, assumes the asymptotic form |^| at the limits of 


integration 
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If the bending stiffness has a negligible magnitude the version of (6) 
for a mass-loaded plate. 


1F„ ~ 




v(x) 2 / A -k e:!p(-lSK) 

2’"”“ - ^ 

/? -k + y 


d? 


‘^0 5, , ^>'"0 ; e.p(-iCx) 

/r-k + y 


(7) 


presents the specific complex-valued function 


I<x) = / d5 (8) 

dealt with by Crighton in the first application of his method for estimating 
Fourier integrals at small values of the parametric variable; and the available 
estimate, which serves for later comparison, is 


I(x) = -21og x-2(log f -^) - 2 ^ sinh“^(^)- iI^-2C +Tryx , x -»• 0 (9) 

n 2 

where y ~ C = 0.5772 ... denotes the Euler constant. 

To arrive at a characterization of the (even) function I(x) suitable 

for development when small values of x are contemplated, assume that x > 0 

and displace the contour into the lower half of the ^-plane where the integrand 

is exponentially small at infinity; then I(x) is expressible in terms of a 

principal value integral along the branch cut between and -k , together 

/2 2 

with a (half-) residue contribution from the pole at ^ = -A +y = -y the 
lower side thereof. Consequent to a change in sign for the integration variable 
along the cut there obtains, accordingly, 


I(x) = 


k ixx. j V “ ixx. j V 
^ e (-dx) ^ ^ e (-dr) 

/T~2 ^ k /TTr ^ 

/t -k + y -/t -k + y 


a- 

Y 


00 
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00 fl 2 

^ r ixx /t “k , . y iyx 

= 2 -f e j- dT - iri e ' 

k T “Y ' 

00 ITX - 00 ixx , 

- 2 S dT + 2U^ / ^ - di Ji e”’‘ . =. > 0 (10) 

/tW /?:? ^ "" 

with the symbol / used to designate a principal value. The first integral 
appearing in (10) can be explicitly described by means of a Hankel function, viz. 


ITX 


/ 


/ 


2 , 2 


dT = ^ (kx) , X > 0 


( 11 ) 


/t -k^ 

and thus the second defines a function 


J(x) = -f 


00 ixx 
e 


dx 


2 2 


( 12 ) 




which satisfies the differential equation 


^ (kx) . 


dx 


(13) 


It is a simple matter to express the general solution of the latter equation, 
namely 


J(x) = Ae^"^^ + Be - ~ f sin y(x-x* )Hq^^ (kx’ ) dx’ 
and the consequent relations 

a 

A+B = J(0) = -f 


(14) 


dx 


Y(A-B) = -iJ’(O) = f 


k / 2 , 2 . 2 2. 

/x -k (x -Y ) 

QQ 

xdx 


(15) 

(16) 


^ v4^-k^ (T^V) 


determine both the constants A,B . Inasmuch as the principal value integral 
(16) vanishes, A = B , and thus 
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A = B = i / - - - 

^ k /T 72 , 2 2, 
✓t -k (t -y ) 




whence 


^ ~ ^ ^ Y(x“x')H^^\kx*) dx' . (18) 

Combining (10), (11), and (18) supplies an alternative representation of 
the function defined in (8), 


I(x) = 77i (kx) - ” log cos Yx - TTi 

0 Y ^ Y-y Y 


Y = v^^+y^ 


- TTi -^ / sin Y(x-x*)HQ^\kx* ) dx' , x > 0 


which lends itself readily to series expansion for small x . The specific 
estimate (9) furnished by Crighton emerges from (19) once the approximations 

„(1) \ • 1 . 2i T kx , 2i p 

Hq (kx) = 1+ — log-j^ + — C 

cos Yx = 1 , e^"^ = 1 + iYx 

2 

are employed and the integral (of order x log x ) is neglected, 

A slightly different procedure for passing between the representations 
(8) and (19) commences with the resolution 

I(x) = / — ^ exp(-i?x) dc = I^(x) - yl 2 (x) + y I^(x) (2 

_oo c - y 


I (x) = f -—nr- dc = 2 / dT = (kx) 


<» -i ?x iyx 

l2(x) = / ^ j dC = -TTi 

-00 C -Y ' 


X > 0 
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and 


I « . 7 „^yp(-isx) — _ 

-eo /TTT .2 2. 

/c -k (? -Y ) 


(23) 


Since 


2 

(-^ + Y^)l3(x) = -TTi HQ^^(kx) 
dx 


and 


laCx) 


= + 3 e - — f sin y(x-x* (kx* ) dx’ 

Y n 0 


(24) 


it follows that 


and 


a +6 = 13 ( 0 ) 


/ 

„JCO 


dZ 


/2T2 ,2 2, 

/? -k (Z -Y ) 


(25) 


Y(a-B) = -il’(O) = - / -- - . - 5 . ^ = 0 , 

^ ll 2 .2 2 . 

A -k (C -Y ) 

the null value of the last integral being a consequence of odd symmetry for 
the integrand. On deforming the contour of the integral (25) around the 
branch cut extending from k to and observing the cancellation of half- 
residues from semi-circular arcs traced in the same sense above and below the 
point C = Y thereupon, the connection 


CO 00 J 

/ is 2 i ^ (26) 

/TT2 /^2 2. k /2 , 2 , 2 2. 

/? -k (C -Y ) -k (t -Y ) 

results; hence (15) and (25) jointly indicate that a = 3 = 2A = 2B and the 
identity of (24) and (19) follows. 

It is appropriate to record, after noting the relations 


/ sin Yx (kx) ^ 
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and 


f cos yx (kx) dx = - Y > k 

0 ^ 


a complementary form of the representation (19) , namely 


I(x) = TTi H^^^kx) + T^± — f sin Y(x-x' (kx' ) dx' (27) 

X 

which permits a systematic approximation for I(x) at large values of x . 

Upon rewriting the response integral (6) , that incorporates both the 
effect of plate mass and stiffness, in the fashion adopted by Crighton, 


v(x) = - 


i^ 

2TTB 


_oo 


(C^-k^)(?^-k^ 

P(?) 


+ yk F(?) 
p 


-i?x j- 
e d^ 


where 


(28) 


P(0 


(?^-k'^)^(?^-k^) - yV 

P P 


5 

n 

n=l 




( 29 ) 


and 


F(C) 


P(Q 


(30) 


attention focusses on the second contribution thereto, since the first can be 


immediately evaluated in terms of the residues at the relevant poles of the 


rational function 1/P(?) • Further to the use of the partial fraction de- 


composition 


1 

P(?) 


5 

Z 

n=l 


2a ^ 
n n 



n 


a 

n 


dC 




n 


5 

Z a ? = 0 

- n n 


(31) 


the function 


G(x) 


f F(C)e"^^^ dC 


00 

f 

—00 


p(?) 


-i?x 

e 


d? 


(32) 
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is linked with previously encountered ones, namely J(x) [cf. 


(23) and (12) ] ; thus 


5 <» JT^ 

G(x) = Z / -“22 — exp(-i?x) d? 


n=l 


n 'n 


-00 ^ _r 


n 


-2 I v(cW) / £ a . <-icx) _ 

„-i "" " - ^ 


(33) 


n 


Suppose that the roots , n = 1-5 , of the polynomial equation P(^) = 0 
lie above the integration contour in the while those of opposite 

sign are located below the contour. Then, if a particular is complex- 

valued, additive and distinct contributions to the pertinent integral in (33) 
arise from a pole (at ^ ) and a branch line (between +oo and +k , 


according as x ^ 0 ), viz. 




f exp(-i^x) dc _ iri _e .or exp(iTx) dx ^ n • 

j 2 2 ^ ^ I ^ ^ J , , o 9 > ^ ^ ^ > 




/t -k ^n 


n 


(34) 


and estimates of the latter function for small values of x can be found in 


a manner described above. The pole contribution is absent if is real and 


larger than k in magnitude, whence 
f exp(-lgx) dC _ 2 / exp(lTx) 


dx 






= 4 cos 5 X ^ 


dx 


X 


cos C X 
n 


n n 


k /2 .2 T -C 0 

/t -k 

log 


on ^ si" C (x-x')H^^^ (kx') dx' 
2 _2 2C „ n 0 


(bi^ /Ci'l jii. 


, X > 0 


X 


.( 1 ) 


/ sin C (x-x’)H^ "^(kx*) dx* . 
^ n U 


(35) 
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2, The Point Source Problem 

A corresponding version of the analysis hitherto given is available for 
dealing with fluid-loaded plate excitation by point forces; in this circum- 
stance the plate velocity resulting from a normal force is 

characterized by the Hankel transform integral 


/Z 2 

icoF 00 /C -k 
v(r) = - ^ f 


(Cr)?dC 


2ttB 


0 .A A. /2~2 .4 

(S -k^)vC -k - yk 


(36) 


P ' ■ P 

where r denotes the radial distance from the point of application of the 
force, Jq identifies the zero-order Bessel function and the other symbols 
retain their earlier significance. Utilizing the same rearrangement of inte- 
grand factors which brings the response integral (6) to the form (28), it now 
follows that 


v(r)=^ Z 

zd - nnnpn U n 
n=l 

iwF oo 

■"ziB / (37) 

where the sum contribution, involving the Hankel function , arises from 

the pole singularities of the function 1/P(C) defined in (31) (namely, at 
the points , n = 1-5 , in the upper half of the ^-plane) and the integral 
term involves the function F(^ defined by (30), (31). The Hankel transform 
of F(^) , namely 


G(r) 


/ F(?)Jn(?^)?dC 

0 ^ 


(38) 


is next rewritten [compare (32)-(33)] as 
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5 ~ ^2 , 2 

G(r) = Z 2a c / "I 2 Jn(Cr)i;dg 

n=l ^ " 0 ° 

^ n 

5 2 2 “ CdC 

= 2 E a^C f —2 (39) 

n=l " " " 0 /TT2 

-k ^ ^n 

whence the point force admittance, or ratio v( 0 )/Fq , acquires the representation 




- E a c (cW) / ■ 1 I . 

n=l " ” ^ 0 /2 2 

/C -k ^ ^n' 


After subdividing the range of integration and employing suitable changes of 
variable it turns out that 


r CdC 1 . CdC 1 . gdg 1 

— — 9 9 ” 9 « T X J . - . 9 9 


^1 ^ f 

r dT . ^ dT 

= J -"-9 9 9 “ 1 / —9 9 9 (41) 

0 T - c + k 0 T + c - 

n n 

and thus the individual integrals of (40) are readily evaluated in accordance 
with the specifications (real, complex) of the ; if is real and greater 

than k , the first (a principal value integral) vanishes and the second equals 


1 .2 i 2 .- 1/2 [ ^ n ”^ 

' n 

A knowledge of the function 


„W- 7^-^ 

» ^ « -'n 


( 42 ) 
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(which has the particular value (41) at r = 0 ) is a prerequisite to the 
eventual determination of v(r) from (37) - (39); and, inasmuch as 


(- 


dr 


+ --^+ C")H(r) = - / 
r dr ^n J 


CO jQ(cr)?d5 


ikr 


n 




( 43 ) 


there follows a representation, 


■r 

TT T It T* ^ 

H(r) = AjQ(C^r) + BNQ(C^r) + f ^ 

-f No(C^r) f dr’ , (44) 

which involves both the zero order Bessel and Neumann functions and is ren- 
dered precise by the assignments 


A = H(0) , B = 0 . 


(45) 


Since finite range integrals occur in (44) the latter permits ready estimates 
for H(r) at small magnitudes of r . 

The point admittance and velocity profile for a fluid loaded plate 
have, as a matter of fact, been determined by Filippi and Saadat (1972), 
independently of Crighton, using a straightforward though lengthier transform/ 
inversion procedure; none of these authors, however, consider distributed 
force plate excitations, whose analysis is taken up next. 


3. Distributed Source Problems 

If the thin plate adjacent to an infinite expanse of fluid, is subject 
to a uniform normal pressure, with the magnitude P^e , over a strip 

domain [x] < jl, -«> < z < its velocity factor possesses the one-dimensional 
Fourier integral representation 
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1 n\P 00 /2 2 

= --5^ / . A . -k sis^ j; , 

™ 4 .4, /TTI .4 ^ 

(C -kp)/? -k -ykp 

which reverts to (6) in the limits Jl ->■ 0 , 2 PqJI -> Fq . On neglecting the 
bending stiffness at first, a simplified version of this integral, 
iP <» /2 T 

/jj) ^ 0 j. /c -k exp(-iCx) sin , 

/ 2~2 ^ 

/C -k + y 

IP 00 lllP 00 

= _0 ^ sin_^ _ Jlo j. e:^(-igc) sin_^ ^ 

S ■’ irnu) /2~T ^ 

/? -k + y 

iPo iyPo 

[sgn(5,+x)+sgna-x)] - _ l(x) 

places in evidence the (even) function of x , 

LO,) . 7 sM;i5s)_ii!^ <,5 

whose counterpart, I(x) , is defined in (8), Pursuant to a resolution of 
the sine function into exponentials, there obtains 


L(x) = -i- / ~ (1-e^^^) exp(-i^x) c 

-CO ^ fl~l 

4 -k + y 

with the component parts 


= Lj^(x) + L2 (x) + L^Cx) (48) 


HU ) ■ i r 

^ 2^-00 ^ pr^ 

/? -k + y 


LzCx) = ^ ^ 

—00 ^ 


/c^-k^ 


C -k + y 


“ j^_^-i?(A-x) 

^3^^^ " ~ 21 E 

^00 ^ 


4^-k^ + U 


✓C -k + y 
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When 0 < X < il all the latter can be inferred from the single function 

dC 


F(a,B) = / 

—00 




- e 


, a, 3 > 0 


( 50 ) 


/22 
/? -k + y 

which, on adopting the transformation scheme of §2, is expressible in the form 


F(a,3) = 2 / 
k 


dx 


Let 


and 


whence 


Write 


/T72 

A -k 


iBx_ iax 2 ‘ 

(.- ) + 2y^ 


dx 


1 i 3 T iax 


. /T — 7 2 2 

'' xw ■' ■ I ' 

+ ,1 X . 

Y 


oo , iax 

P(a) = / ^ ^ 


dS - - 2 Hq ’ 


Ct / T \ IcOt /"IN 

P(a) = P(0) - ^ / H^^^ka')da' = ^ [1 - f d?] 

z 0 ^ 0 


dT 


1 e 


iax 


Q(ot) — ^ _ ■ 9 2 T 


and (cf. (12), (18)) 


(51) 


(52) 


(53) 


(54) 


= - Tpi— log cos yOL + 7^ -fsinYOt f cos Yx'H^^^(kx') dx' 
da 2 yp ^ y-\i 2y ^ ' q 'O'' 

^ ( 1 ) 

- cos ya /sin Y^'H' (kx') dx 

0 


•} 


SO that 


j /I \ 

i ’ 2tH .. ■ /•«(!), 


Q(a) = Q(0) - log ^ sin Ya + S H' '(O d? 

2y y 2kY 0 


IcCt 

— f cos Y(a- |)H^^^(0 dC . (55) 

2kY 0 
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On the basis of (51) - (55) a finite integral characterization of F(a,B) 
is reached, viz, 

h6 

F(a,6) = - ^ / (Od? +TTi “ i 1°8 (sin yB-sin ya) 


y ka 


2 

-TT / COS y(3-f)H^^\o d?- 7 cos y(a-|)a5^\?) d?] 

ky 0 ^ ^ 


k^^O 




and thence follows for the bounded velocity v(0) along the mid-line of the 
strip 


v(0) = 


^ _ IUPqL(O) _ ^ ^ 

mw TTmoo mO) TTmOi 


F(0,JI) 


^^0 ^^0^ (1) ^0 iy5/ ^0 Y+U 

= — + — ^ / h;-^^(^) d? ^ (1-e^n + ^ log ^ sin yZ 

hkjo 2 / 0 ^ ^ 2 ^ 2 ^ y-\i ' 

mo3Y 0 iiujOY rnncjoy 

^ ^0 £ (1) 

+ ^ f cos Y(^l-f)H^ dC • (57) 

mky GO 0 

It is also possible to express, in analogous fashion, the mean velocity over 
the width of the strip, 

1 ^ ^ 

^ / v(x) dx -T— 5 — / L(x) dx 

2Z moo 2TTilmG0 


+ y 


1 

xyPo 

SL 

f 

-a 

moo 

2TTilmG0 

l!o_ 

iyPo 

% 

f (■ 
-a 

mco 

TT^moo 


having regard for the result 


s i il.. ;^ )2 dg = k I ^ y jj(l) (2k«,)| + TT-^ 

^ ,/T^ y^ I 0 ° J 




-i ^ log ( 1 -COS 2y£) 

2y 2y-^ 

f sin y(2£ - d? 

2ky-^ 0 ^ 


+ iTT 
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and this enables a characterization of the distributed force admittance/ 
impedance of the plate. 

For points on the plate external to the strip, say x > 5, , the relation 


L(x) = F(x-5,,x+Jl) 


(60) 


obtains and, applying (56), this yields 

L(») . (O dC +-% 

2y k(x-£) ^ 2 y 


log (sin y( 35+^') - sin y(x-Z)) 

2y 


(61) 


2 k(x+il) P k(x— S,) p y. X 

+ iTT { / COS Y(jH-il-|)H^-^^(Od? - / cos Y(x-ii-^)H^ ^(Od?} 


2kY 0 


0 


In the limit 5, -»■ 0 , 

L(x) ^ H-l L)^^o + 0(il2) 

= il{i ^ (kx) - i ^ ^ log ^ cos Yx + i ^ (kx) 

-1-^ / sin y(x-^)Hq^^( 5) d?} 

= ill(x) + 0(5.^) 


and the recovery of the line force expressions (7), (8) from (46) is immediate. 

When the applied pressure acts on a circular area (with radius = a) the 
velocity distribution corresponding to (46) takes the form 
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v(r) = 




m03 




c 0 


J_(Cr)CdC 


f J„(5DJ^(;a) dc - ^ / -±=J d; 


« 


(62) 


fx, -k + y 

as a function of the radial distance from the center of this domain. The 
simplification of (62) commences by noting that 


/ jQ(cr)J^(Ca) d^ = l/2a , 0 < a = r 


0 , 0 < a < r 

l/2a , 0 < a 
1/a , a > r 


(63) 


and adopting the rearrangement 

~ jQ(;r)J. (?a) ~ jQ(Cr)Jj^(Ca) ^ ~ jQ(?r) Jj^(?a) ~ jQ(^r)Jj^(Ca) 

/ — dc = / nzn +vi / — dc - y / — o — 0 


0 /2 , 2 . 

/? -k + y 




0 /2 , 2,2 2. 
Vr, -k (C -Y ) 


n .2 2 

0 5 -Y 


(64) 


for the second term. Let 


~ J»(Cr)J. (Ca) ^ a 2ir 

M(r) = / ; dC = -5— / r'dr' / di|i ® 

0 ° 0 


and define 


xp(ik^r^+r'^-2rr*cos ijj) 
/2 2 

/r +r' -2rr* cos ip 


(65) 


whence the relation 


00 Jn(;r)J (?a) 

N(r) = / -y ± dC 


( 66 ) 


0 Z2 , 2 ,2 2. 

/? -k (C -Y ) 


(-^ + i ^ + Y^)N(r) = -M(r) 
dr 


follows, with a consequent representation 
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r r 

N(r) = AJpCyr) + j J^Cyr) / r'N^ (yr ’ )M(r' ) dr' -j N^Cyr) / r' jQ(yr' )M(r' )dr' 


wherein 


A = N(0) = f 




M. 


2 2 * 




Since 


• d^. 1 d , 2 1 ,, r 

da a 


0 


, ^ ± , ika , V 

dc = (a -1) 


it turns out that 


(67) 

(68) 


A = oJ^(ya) “ '^iWa) / N^(yx) (e^^”^-!) dT+-|^N^(ya) f J^(yx) (e^^"^-!) dx 


= oJj^(ya) 1- (1-e^^^) + ~ J^(ya) f NQ(yx)e^^"^ dx 


kay“ 0 

^ HT /. \ /• T ikX 


- N^(ya) f dT 


0 


(69) 


after invoking some Bessel function properties. A vanishes, independently 
of a , when a = 0 and 

1 = A 


n, JL = = - r ^ 

“ 2 Ma''a=0 2 J nr—r. Jl 2 


0 ^2 5-Y 


SO that 


1 r r dx ^ ^ dx 1 i ^ -1 k 

“ "7 TT- ^ / X^’ ■ ■ w ' d 

^ 0 T -y 0 T +y ' 

and the determination of A is complete. 

The third and last integral in (64), say 

00 J (?r)J (ca) 

P(r) = / - ^ 


0 


2 2 
C -Y 


d? , 


can also be recast by means of its differential equation 


(70) 


(71) 


(■^ + 7 ^ + Y^)I’(r) = - / JQ(?r)J^(?a) d? = - 7 , a > r 
dr 0 


( 72 ) 
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and this yields 

r r 

P(r) = PJ^Cyr) + '^ ^ r'NQ(yr') dr' - ~ N^Cyr) f r'J^ (yr') dr' 

= 3JQ(yr) - 

ay 

or 


where 


P(r) = P(0)jQ(yr) - (l-J^Cyr)) 

ay 

°o J (Ca) d? oo 

P(0) = / i / J,(ya) d, 


(73) 


(74) 


0 ^ -Y '0 

and a principal value is understood. Differentiation with respect to a 
in the formulas (cf. Luke (1962)) 


{ in - -f [H„(Ya)-N„(ya)l 


and 


CO J^(ca) d? ^ 


/ - = t [H^(ya)d-N^(ya)] , 


where Hq denotes the Struve function, yields 


CO J (^a) 
/ ^ 


0 




and 


CO J (ca) 

^ -^dC=:^ + f [H'(ya)+N;(ya)l 


0 --Y 

(with primes indicating an argument derivative) , whence 


(75) 


(76) 


P(r) N^(ya)jQ(yr) , r < a . (77) 

ay ^ 

Collectively, (65), (67), (69), (70), (71), and (77) provide an efficient 
basis for estimating the last integral in (64) — and thus the normal velocity 
v(r) — within the area directly subjected to external surface pressure; and 
they furnish, in particular, an expression 



- 20 - 


v(0) = 


iPp ^ layP g ^ ^ ^ 2 

mo mu) 

IF, 


+ y A + -^ + ^ N^(ya)} 


ay 


“ ^ ^ (e^^^-1) + i ^ tan"\^)J (ya) - ^ ^ N^(ya) (78) 

y I M Y 2y 

- ^ 7 "o'k 1?^ ’'i‘^"> 7 •'o'J 

for the velocity at the center thereof. The mean value of v(r) on this area 


V = —r f rv(r) dr 
a 0 

TOO mto C 


dC 




(79) 


-k + y 

can be related to previous integrals whose display in convenient forms is 


available; and the link appears after a sequence of manipulations 


“ J, (?a) 
/ 


dC 


, ^/2 ~ J (2?a sin 0) , 

= 4 / de / ^ 


/e-i^ 


-k + y 
TT/2 CO 


9 “ f 1 ,1 ‘Jl 

■?/ 


(2^a sin 0) 


-) 


1 




dC 


, 7t/2 

= 4 / d0 
^ 0 


2(y-ik) 2a sin 0 


^ <=o j (2 sin 0) 

f ^ 


dC 


0 ^ 




1 1 d 6 r ®> cVW-2u^ .. 

*7 r\ J TOO 


2 y-ik ira p"' sin 0 p"' ^ — 2 '-2 2^2 

/c -li 


1 1 _ d0 

2 y-ik *^^0 sin 6 


r-k^ 

00 J (2?a sin 0) 


u -y y 


I r-1 

i n 


-dL 


0 AW 


^2 2 

? -y 




00 (2^a sin 0) 


y dy 




dC _ Jii 

W ^ 0 


00 (2?a sin 0) 

/ - 


2 2 

c -y 


d? 


( 80 ) 


which turns up representatives of the integrals defined in (68) and (71) . 
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There is, likewise, sufficient information at hand to obtain forthwith 
the corresponding measures of the plate velocity when bending stiffness 
plays a role. The integrals 


7 no dt - 2 L c (cW) 7 

^ C T n n n ^ 

~oo n=l -0® 


and 


( 5 ^- 4 ) 


“ Jd (Ca) 5 _ „ 00 j (^r)J (?a) 

f ^ — J (Cr)?dC = 2 Z a c (C-k^) / 

0 ^ ° n=l ^ ^ ^ 


0 2 ,_2 _2 


r-k" (? -?;) 

in particular, make an appearance and the above presentations reveal that 
functions of the type (54) and (66) are involved. 
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